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Introduction 



Quasi-one-dimensional (quasi-lD) quantum gases can be formed in prolongated potentials (atom waveguides) when- 
ever the atom energies are much less then the transverse excitation energies and only the ground transverse state is 
t~^. , substantially populated. Such gases have been recently realized in 2D optical lattices pi 0, 0, 0, 0, 0, 0, 0] , elongated 
(N 1 atomic traps %M, EI El El 111 , and atomic chips 0, El E3, El, El • 

Properties of quasi-lD gases are drastically different from the 3D ones (see review [2(| and the references therein) 
Under certain conditions quasi-lD gases can be described by the integrable Lieb-Liniger-McGuire (LLMG) model 
O ^ [2lL l22l. [23!] of ID indistinguishable Bose atoms with zero-range energy-independent interactions. This model has an 
exact Bethe-ansatz solution and can describe di- and multi-atomic molecules in the case of attractive interactions. 
However, reflection and dissociation in atom-molecule collisions and three-atomic association are forbidden within 
the LLMG model. Its high internal symmetry can be lifted in more realistic systems, leading to observable effects of 
integrability, such as atom-diatom reflection for non- identical atoms [11 . Richer physical phenomena can be provided 
by resonant interactions. 

• Resonances in interactions of quasi-lD atoms can arise both from virtual excitation of transverse modes 
(confinement-induced resonance, or CIR, [HHH) and from atomic internal structure (Feshbach resonance, see [27l| ) . 
Combined effects of Feshbach resonances and CIR on two-atomic systems were analyzed in [H , where parameters 
of ID scattering were related to parameters of 3D scattering and harmonic waveguide. Effects of an anharmonic 
Q_i, confinement, finite-range interactions, and other additional effects were analyzed in (30l. l3lL l32l. [33l l3~il l35l [3r| . 

Resonant interactions can lift integrability already in three-body systems. CIR can lead to atom-diatom reflection 



37 1 (a similar process of atom-soliton reflection was analyzed in [381 ]) and to thermalization in three- atom collisions 
39(. Narrow Feshbach resonances can pro vide richer physical phenomena, such as association in three-atom collisions, 
dissociation in atom-diatom collisions [4d| | and stabilization of Feshbach molecules [4l| . 

An effective approach to three-body problems was developed on the base of Faddeev equations (see books [H|, HU 
and the references therein). Exact analytical [3, El and numerical [46| solutions of Faddeev equations for the LLMG 
model agree with the Bethe-ansatz solution. Numerical solution of Faddeev-Lovelace equations has been applied in 
C ' [11 to asymmetric systems and in [47[ to ID systems with generic separable interactions. An equivalent method was 
used in [37| . An alternative method of numerical solution of the Schrodinger equation in hyperspherical coordinates 
was applied to ID problems in [48]. Three-body bound states for resonant interactions were analyzed in [4^ ] using a 
variational approach. 

The present work gives a comprehensive description of the Faddeev approach to ID problems, generalizing the 
k^j ' approach of [4(il . |4~0 ] to multichannel resonant interactions. The second-quantized Hamiltonian for this case is in- 
troduced in Sec. II Al Two-body scattering and bound states are analyzed in Sec. II Bl Section II CI describes the 
three-body multichannel problem. Asymptotic states and transition amplitudes for this problem are introduced in 
Sec. II DI Section ITEl relates the multichannel transition amplitudes to solutions of the Faddeev-Lovelace equation for 
an effective single-channel problem. Probabilities of bound-bound and bound-free transitions are presented in Sec. 
II Fl Section II Gl describes special cases of the two-channel system and deactivation problem. The solution method of 
Faddeev-Lovelace equations is described in Sec. [Ill The present regularization procedure is applicable to generic ID 
problems with separable interactions. The singularities related to the poles of two-body T matrix, which are similar 
to the 3D case, are treated here using a generalization of the approach [11. The present approach also regularizes 
the singularities of the free-atom Green function, which are different from the 3D case. They appear in the processes 
of dissociation and association and were not considered in [24j . Derivation of certain properties of resonant two- and 
three-body models is included in the Appendices. 
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I. ONE-DIMENSIONAL BOSONS WITH MULTICHANNEL ZERO-RANGE INTERACTIONS 



Hamiltonian 



Consider a gas of ID Bose atoms described by the annihilation operators if? a (z). The model includes several two- 
body channels described by the annihilation operators \fr m (z) of molecules with the energies D m . The Hamiltonian 
of the system has the form (using units with Plank's constant Ti = 1), 



H 



dz\¥ a (z) 



2m dz* + 2 *«W W «W 



*»(*) 



4m <9z 2 



An] *m (2) + E [ft™ ( Z ) + ( Z )] + E I' ' 



(i) 



where m is the atomic mass and U a is the non-resonant interatomic interaction strength. At negative values of U a 
the open channel, formed by two atoms in the initial state has a bound state too. The interactions 



Km (z) = g m ¥ m (z) # a (z) * Q (z) 
couple the open channel to other (molecular) channels and 



V mm > (z) = dmm'^m (z) * m / (z) . 



dm'm — d mm ,, 



(2) 



(3) 



couple the molecular channels [the matrix d m / m will be completed by diagonal elements, see Eq. © below]. The 
system with a single molecular state was analyzed in [28l . |29| , where the one-dimensional parameters were related to 
the atomic collision and waveguide parameters. Additional molecular states with large negative D m were introduced 
in ED. 



B. Two-body system 

A state vector of the two-atom system can be represented as a superposition of atomic and molecular states, 



00 

* 2 > = / dz c e iP *< [i= J dz V W (z) *t (, c _ f) $t ( Zc + |) + £ <pl<»W m (z c ) 



I vac), 



(4) 



where z c and P are the center-of-mass position and momentum, respectively, z is the interatomic distance, and |vac) 
is the physical vacuum state. Substitution of the state vector ^ into the Schrodinger equation with the Hamiltonian 
([1]) leads to the coupled equations for the channel functions ip^ (z) and ip( m \ 



EipW (z) = 



1 d 2 

+ 

m az z 



,(0) 



{z) + V28(z)Y J 9l. 



(m) 



Eip^ = D m9 ^ + V2. 9 „^ (0) (0) + J2 d-mm'^ 

Here E is the energy in the center-of-mass system. 
The molecular channel functions can be expressed as 

= yfapM (0) G m (E) , 

where the vector G m is a solution of the system of linear equations 

2_j drum' (E) G m > (E) = — g m 

m 1 

and the diagonal elements of the matrix d mm i (E) are defined as 

d mm (E) = D m -E-iO. 



(5) 
(6) 



(7) 



(8) 



(9) 



3 



If g m — —S mm i , then G m (E) has a physical sense of a Green function, and the infinitesimal imaginary part of d mm (E) 
specifies its retarding behavior. 

Elimination of the molecular channels leads to a single equation for the open-channel function, 



--^ + U eS (E)8(z) 
m az z 



EipW (z) = 

where the effective energy- dependent interaction strength 

U eS (E) = U a + 2j29* m G m (E) 



(10) 



(11) 



incorporates effects of all channels. In the case of a single molecular channel it is reduced to the expression for U c s (E) 
in Hi]. 

Obviously, the wavefunction of a two-body scattering state can be expressed as 



<pM (z) = (2^- i/2 



ikz ■ rji j ] \ ik\z\ 

3 -t—T lD (k)e 1 1 



in terms of the two-body T matrix 



Tid (k) 



Uj (k 2 /m) + -mk- 1 



(12) 



(13) 



which depends on the relative momentum k of two colliding atoms (E — k 2 /m). Tid (k) has poles on the positive 
imaginary axis, k = iK n , where n n are solutions of the equation 



'-U eS 



K 

m 



(14) 



The poles correspond to two-body bound states (diatoms) with the energies E n = —K^/m. The diatom states are 
orthogonal (see App. [X| and the orthonormality conditions have the form 



(15) 



Solutions of Eq. {5| can be expressed in terms of det(d), the determinant of the matrix d mm i (k 2 /to), which is 

a polynomial of the degree 2M mo i in k, and the cofactors ddet(d) /dd m m', which are polynomials of the degree 
2 (M mo i — 1) in k (where M mo \ is the number of the molecular channels). Using Eq. (jTT|) . Eq. IT4]) can be represented 
then in the polynomial form, 



— K + U a ) det(d) - V 
m J z — ' 



9m' 



m,m' 



ddet (d) 
dd„ 



-071 



0. 



(16) 



This equation has 2Af mol + 1 roots (only the real positive roots correspond to the bound states). Thus the energy- 
dependent interaction can lead to multiple bound states in spite of its zero range. The diatoms are superpositions of 
the open and molecular channels. The open-channel component can be expressed as 



(z) = ^ (0)exp(- K „|z| 



,(0) 



(17) 



where ipi^\o) is determined by the normalization conditions (|15p . and the open-channel contribution to the bound 
state is 



dz\<p™ (z)\ 



l + 2j2\G m (E n )\ 2 K r , 



(18) 



Although the present model approximates the molecular-channel states to be infinitesimal in size, the diatoms have 
finite sizes (~ wC 1 ). 
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C. Three-body systems 



A system of three atoms is described by the state vector 



l* 3 > = 



)= [ d 3 z^ (z l7 z 2 , z 3 ) tft (zj) *J ( Z2 ) $t (Z3 ) + j dzdz m ^ (z, z m ) ¥ a (z) & m (z m ) 



V6 



|vac), (19) 



leading to the coupled equations for the wavefunctions of the three-atom tfj^ (z%, Z2, 23) and atom-molecule 
-0 ( - m ^ (z,z m ) channels. Their momentum representations, ffi ' (q,k) and (q), are defined by 



V> (0) (zi,z 2 , z 3 ) = (2?r) 3/2 J dqjdkjip^ (qj, kj) exp (iqj (zj 



ikj (zj> - Zj»)+iP 



Z1+Z2 + z 3 



(20) 



^ (ro) (z, z m ) = (27T)- 1 y exp (i Qj (z - z m ) + iP 



.„z + 2z r< 



where P is center-of-mass momentum and qj, kj are momenta in any of the three sets of Jacobi coordinates (j = 1, 2, 3). 
The momentum qj is the relative momentum of the j th atom and the center-of-mass of the j' th and j" th atoms, 
while kj is the relative momentum of the f th and j" th atoms, where j, j' , and j" are cyclic permutations of 1, 2, 
and 3. The momenta in different Jacobi coordinate sets are related as 



1 



1 



(21) 



These relations correspond to permutations of the atoms. In what follows the center-of-mass system (P = 0) is used. 
The channel wavefunctions in the momentum representation satisfy the following equations 



3 
4m 



q 2 + ±- k A^ (g . )fe . ) + ^ f ^(o) {quk[) + i:i « ) -^Y J 9: n ^ {m) i<i l ) 

' l=\ J m 1=1 



(22) 



4m 



Urn) 



1/2 



g m j dkj^ 0) (qj,kj)+ 4.m'^ m '' (ft) 



where E is the energy in the center-of-mass system. These equations are invariant over transformations (|21[) and can 
describe indistinguishable Bose atoms. 

Equations (|22[) can be represented in a compact vector form. Consider for the moment distinguishable atoms and 
introduce the vector wavefunction tp — {tp^ 01 '}, where a has values 0, ml, m2, and m3 for all m. The component 
^p(mj) — ^>( m ) (g^) describes the atom-molecule channel with free j th atom. Equations (|22|) can be then written out 
in the form 



where elements of the matrices Hq and Ui are defined, respectively, as 



(23) 



j_ra a 

H Q - 



(«) 



1/2 



^ Qm $ct' ,mj $lj 



(37r) 1 ^ 2 5 , m ( ^a'o+ £ d m i m 6 a i <m 'j'6ij' 

m',j' 



S a0 J dk$W 



(24) 



(25) 
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The definition of H$ a is unambiguous since the kinetic energy -g^qj + -^kj keeps the same value in the three Jacobi 
coordinate sets. 

The atom-molecule channel wavefunctions can be expressed in terms of the three-atom one, 

(g .) = Q V2 Gm (e l_ q f\ J dk .^o) {qhkj) > (26) 

where G m is defined by Eq. ||5J). Elimination of the atom- molecule channels leads to a single equation for the three- 
atom wavefunction, 

fe, %) = (^ 2 + ife]) ^ („,*,) + £ L>r , (27) 

where the interaction operators are defined as 



1=1 



Uf 



= ±-U eS (e - J-g? ) / dk[^ fa, k[) (28) 



27r V 4m 



with J7 c ff given by Eq . fllip . Elimination of channels and energy-dependent potentials in 3D three-body problems were 
considered in [50, [5l|, [52j . 

However, in three-body problems effects of the eliminated channels extend beyond the interaction strength. The 
additional effects are analyzed below. 

D. Multichannel Faddeev approach 

Equation (|23p describes a vector three-body problem. Similar scalar problems were extensively studied using 
Faddeev equations (see books [42l |43]|). In line with this approach let us introduce vector functions Xoq and Xinp, 
describing the asymptotic channels. The wavefunction of three free atoms with the momenta Q%, Q2, and Q3, such 
that Qi + Q2 + Q3 = P = 0, has the form 

*S fe, fcj) = S a0 5 {q 3 - S (ki - Qy 2 9j " ) (29) 

in each of the Jacobi coordinate sets. The asymptotic function containing the free I th atom and a diatom in the state 
n is 

xZ (©,**) = « (« - P) ft) , x£° (®) - M (« - P) ^ m) ; (30) 
where p is the relative momentum of the atom and diatom and 

is the momentum representation of the two-body wavefunction (|17[). The asymptotic functions satisfy the following 
equations 

-E-Xoq = -ffoXoQ 

(32) 

-^Xinp = ( HO + Ul ) Xinp 



where E = Q 2 / (2m) in the first equation and E — 3p 2 / (4m) — K n /m in the second one. 
The elements of the scattering matrix can be expressed as 



Si'n'p'inp — S [p' —p)6 nn >5u> — 2ni8 3 — - Xi> n > p i i np 

\ 4m m / 

(33) 



Q2 p 2 K 2- 
SoQ.lnp = —2-KiS 3- ! XQQjnp 

1 2m 4m m ' 
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in terms of the transition amplitudes 

Xl'n'p',lnp = (Xl'n'p' \ £ ^"l^inp) (34) 



3 



XoQMv = (Xoq\ ^i"l^/np> (35) 



l"=l 



Here the scattering-state solutions i/ji np of Eq. (|23|) have xinp as the incident wave. 
Given energy E and internal state n, the momentum p can have values ±p„, where 



Pn = 2^(mE + K*) /3. (36) 
The bound-bound 5 matrix elements can be then represented as 

4l7TZTTI 

Sl> n >p>,ln±p n = 8 (V TPn) <W<*H' 5 [<$ (p' - Pn') + 5 (V + Pra')] X l'n'p' ,ln±p n , (37) 

"Pn' 

describing transmission and reflection in atom-diatom collisions with possible rearrangement of atoms and change of 
the diatom internal state. 

Indistinguishable Bose atoms are described by the wavefunctions ip^ (^1,^2,23) which are independent over per- 
mutations of each pair of the atomic coordinates z\, Z2, and 2:3. Then the momentum-representation wavefunctions 
ip( > (qj, kj), ^A m ) (c[j) are invariant over transformations ([21) . and the components are independent of j. A sym- 

metric scattering-state solution xp np = ^ ipinp is obtained with a symmetrized incident wave 3 -1 / 2 ^2 Xinp since 

1 1 
Eqs. (|22p and (|23[) are invariant over transformations (|2ip . Then a symmetric atom-diatom collision wavefunction in 
the coordinate representation [see Eq. (|2"0|) ] has the asymptotic 

o-i/2V~^ Airim 1 , , 

tynpn ^ 3 / Xlnp n ~ / ( 3 n 'Pn' M l Pn'X.jn' Pn i T" Xjn'—p n i,ln Pn Xjn'— Pn i ) 

« , Pn' 



£=1 



-27ri y d 3 Q<5 (Qx + Q 2 + Q 3 ) <S - X q,;„ p „x q 



(38) 



where the three terms in the square brackets correspond to the incident wave, transmission-reflection, and dissociation 
of the diatom, respectively, and x are the- coordinate representations of the functions X- 



E. Effective single-channel problem 

Taking into account Eqs. (|25|) and pO]) , and expressing ^j"™' 7 '' in terms of ip^ using Eq. ([26]), one can represent 
the transition amplitudes in the form of 



X lw±Pn , Mpn = (x$L Pn , I £ tWgJj ( 39 ) 

*0Q,*np„ = <*{jg| E ( 4 °) 



;"=i 

involving only the three-atom components of the wavefunctions. These components are proportional to solutions 

of the effective single-channel three-body problem with interaction J7f ff . Indeed, consider the solution of the scalar 
equation 

Exm P (qi, h) = H^Xlnv («, h) + (if [xm P ] , (41) 

of the form 

) fqf + fcf-r^ ' (42) 
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where only the normalization factor depends on n for fixed p. The solutions are normalized as (xinp'lxinp) — 
AS (p' — p), where A = 1 for p = p n . Taking into account Eqs. ([3"0[) . ([3"Tj) , and l|56| ). the wavefunctions in Eq. 
(|39| can be related to the single-channel solutions for an arbitrary diatom state n = 0, chosen as the reference state, 



(43) 



Here V'iOpn is the scattering-state solution of Eq. (|27|) with the incident wave xzop„ (Qu h) which is off-shell [Sp^/ (4m) + 
Kq/tti ^ E for n^O]. Therefore, all transition amplitudes in Eq. |38|) can be related to off-shell amplitudes for the 
effective single-channel interaction Uf and a single diatom state n = as, 

Xi> n >± Pn ,,i nPn = \/W n ,w J Kn ' K f (Xl'0± Pn ,\ E Ui?\'4>l0p n ) 

(44) 

X 0QMp = (^fi) (xoq| E ^rl^opj 
V K " / ;»=i 

[Xpq = xqq since Xoq nas the three-atom component only, see Eq. ([25]) ]. Thus, solutions of the effective single-channel 
problem (|27j) provide all necessary information on multichannel scattering. 

describes a three-body problem with separable two-body interactions. Following Lovelace (see 
43j and an adaptation of this approach to ID problems in [47]). the problem can be reduced to an 



Equation ([27 
also books [4j 
integral equation 

ui^ f 

X (q 1 , q ) = 2Z (q\ qo) + 7TH d 1 Z «) T io (* (9)) * 9o) (45) 
for the symmetric transition amplitude 

3 

*(</,<?) =E<Wl E Wn g ) (46) 

(it is independent of I', see [42L I43I I53I] ) . Here 

Z V, 9) = — F ^- n 1 * ( 47 ) 

7rrn mi/ + il) — q z — qq — q' z 

is the Green function for three free atoms, expressed in terms of the momenta q in two Jacobi coordinate sets, the 
two-body T matrix T\r> (k) is given by Eq. (|13[) . and the relative momentum of two atoms is expressed in terms of q 
as 



k (q) = ^mE + iQ-ZqZ/A. (48) 



F. Transition probabilities 



Identifying the transmission and reflection amplitudes in Eq. ([38[) and using Eqs. (|44[) and ([46)) one can express the 
reflection and transmission probabilities in terms of the symmetric transition amplitude 



f _ 16tt 2 m 2 W n ,W n a 

r n>n ~ a 6 V Pn' , Pn 



(49) 



47rm 
1 _ i Wn 

3Pn 



X (Pn,Pn) 
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Atom-diatom collisions can also result in dissociation of the diatom. This process is described by the last term in Eq. 
(155)) . However, the total dissociation probability is related to probabilities of other processes as 

^ diss = i-E( p «'« n + p « ? «)- ( 5 °) 

n' 

Diatoms can be formed in three- atom collisions. The probability of association per unit time can be found, like in the 
3D case [13] , from the equation 

2 2 

°L _ 
4m m 



C s ° c (Q) = 2vr| ^|<Q| E Ui\n,p)\ 2 S (s^-^-E). (51) 



Here 



iQHffi-^xonCQ) (52) 
n 

is the wavefunction of three free atoms with momenta Qi, Q2, and Q3 (Qi + Q2 + Q3 = 0), and the energy is 
E = Q 2 / (2m). The wavefunction (|52[) is normalized per unit in the box of the length L and symmetrized over all 6 
possible permutations II (Q) of the momenta Qi. The wavefunction of the interacting atom and diatom in the state 
n, 

\n,p) = ( wf I E^«p> (53) 



1=1 



has the box normalization too and is symmetrized over the atom permutations. 

In the system of N a t atoms with momenta Qi the rate of formation of molecules in the state n can be evaluated by 
summation over all trios 

dN n _ \ -* assoc / n Qii + Qi 2 + Qi 3 n Qji + Qi 2 + Qi 3 n Qix + Qii + Qh \ / KA \ 

~dT ~ 2^ w ™ 3 3 '^3 3 ]• (54) 

Replacing the summation by integration with the atomic momentum distribution / at (Q), one obtains the rate equation 
for the molecular density p n = N n /L, 

d Pn a /" ^nu fn Q1 + Q2 + Q3 n Q1 + Q2 + Q3 n Qi + Q2 ±Qs \ f (n w rn w rn x 

-fa = Pat /« V^3« \ Q\ 2 3 3 J /at (Qi) /at (C^2) /at ((^3) , 

(55) 

where /9 at = N ai /L is the atomic density and 

#3« (Qi, Q2, Q 3 ) - ^< soc (Qi, Qa, Os) 



2--/// :l /,-' ! 



X] ^ ( fc (Qi)) * (Q;,p») I 2 + 1 E T ^ ( fc (Oi)) x (Qi. -p») 

i=i i=i 



(56) 



is the association rate coefficient. [Here p n is given by Eq. ([56"]) and Qi + Q2 + Q3 = again.] This expression is 
derived using Eq s. (f5"Tl) , ([52^1 . (|55|) . (|55|) . (|31)). and the general relation between bound- free and bound-bound transition 
amplitudes (see |42|. |43I|). 

3 3 

(XOQl E W«*J = -y (2^^)" 1/2 E ^ ( fc (Qi)) (XiOQ, I E W«*J ( 57 ) 
/"=i i=i z"#i 

The rate coefficient is defined according to conventional chemical notation and is three times larger then the K3 
used in 
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G. Special cases 

1. Two-channel case 

In this case (4p| . when only one molecular channel (to = c) is present, the vector G m (E) contains only one element, 



GAE) = p . (58) 

y ' E + iO-D c y ' 



Therefore, the effective interaction strength is 



U eS (E)=U a + F ^ c|2 „ (59) 
E + tO — D c 

and Eq. (fTo) for the poles of T 1D is reduced to the cubic equation [H, |49| 

k 3 + ~^U a K 2 + mD c K + ^m 2 D c U a - m 2 \g c \ 2 = 0. (60) 

It has two positive roots Ko.i whenever U a < and D c < 2\g c \ 2 /U a , and one positive root k otherwise. The 
open-channel contributions to the corresponding bound states are [29( , 

Wn = <+ mD c (61) 

" 3k,2 + mU a K n + mD c 



2. Deactivation case 



The general problem with multiple molecular channels can be substantially simplified in the case when one of the 
molecular channels, the closed channel (to = c), is coupled to all other channels, while other molecular channels, the 
deactivation products (to 6 {d}), are not mutually coupled (d d d' = 0, d ^ d'). In this case the effective interaction 
strength can be expressed as 

U cS (E) = U a + 2Y ^ + 2(\g c \ 2 + g* c y dcd9d )( E + iO-D c -Y ^ ) . (62) 

V ' ^E + iO-Dd y ycl yc ^E + iO-DdJ\ ^E + iO-Dd) V ; 

If the deactivation product states lie far below the open channel threshold, i. e. Dd < and 

\D d \ ^max(\d 2 cd /D c \,\g d d c d/g c \,\9 2 d/Ua\,\E\,\Dc\,7nU 2 ), (63) 

the effective interaction strength is approximately the same as in the two-channel case (|59|) . and Eq. ([14]) has roots 
which are approximately determined by the same cubic equation (|60j) . 
Equation (TT4")) also has roots 

K d « y/m\D d \, (64) 
one per each deactivation product state. The open-channel contributions to the corresponding bound states are 

W d tt\^\g d \ 2 \D d \-*' 2 . (65) 

The rate coefficient for deactivation of the diatom state n onto all states {d} can be expressed, using Eqs. P9"| . (163)) 
and (|C3p . as 

K ^ = ^ E W + PdD = E (°- °' °) + 7^1 (0, 0) | 2 , (66) 

where the coefficients ^ ad and 7 C( ; are expressed as 

/3 5 4to\ 1/4 /3 3 m\ 1/4 
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and the symmetric wavefunctions 



CI (0,0,0) 



2^ 



im 2 K n 
4«o 



dq 



Tid (k(q))X (q,p r , 



41(o,o) = -^ 



W„K T . 



toZX 



+ 



2kq 



k(q) 

Tip (k(g))X (q,p n ) 
U a (fc 2 (q) - mD c ) + 2to| 5c | 2 



(68) 



are calculated with the symmetric transition amplitudes X (q,p n ) and the two-body T matrix T\u (k) for the two- 
channel case (see App. [Cjl . 



II. REGULARIZATION OF ONE-DIMENSIONAL FADDEEV-LOVELACE EQUATIONS 



The previous section demonstrates that a solution of Eq. (|45[) provides all necessary information for evaluation of 
transition probabilities. This equation is a particular case of the general Lovelace equation (see [42|, HH, |U [Hj]) for 
separable potentials Vv* (q) v (q') with the interaction strength V and formfactors v (q'). For zero-range potentials 
the formfactors are constant functions, v (q) = (2n) , in the momentum representation and v (z) — S (z) in the 
coordinate representation [indeed, Vv* (z) J dz'v (z r ) <p (z r ) — VS (z) tp (z) in this case]. In the general case the two- 
body T matrix is expressed as v*T\dv in terms of the function T±d in Eq. (j4"5|) and the free Green function Z (q', q) 
has the form 



*tf.d-g y ' , r^ ( -ff-'? . (69) 

mh + lU — q — qq — q 



As in the 3D case (see [i2|), solution of ID Lovelace equations can be troubled by singularities. The regularization 
procedure described below is applicable to the general case of separable potentials. 

The singularities can arise from both the two-body T matrix and the free Green function. Poles of T\d (k (q)) 
at q — ±p n are related to two-body bound states and appear for each scattering process, whenever the energy lies 
above the lowest bound state, E > — max (k 2 ) /to [see Eq. (fBl)])]. The free Green function Z(q',q) has singularities 
at q = q± (q') = —q'/2 ± k (q / ), where the function k (q) is given by Eq. (|4"8]) . Since q' and k (q') are the momenta q 
and k in one of the Jacobi coordinate systems, Eq. (|21[) allows to recognize q± (q') as the momenta q in other systems, 
where the atoms are permuted. Therefore the poles correspond to non-diffractive collisions where the momenta keep 
their values and are only exchanged by the atoms. Then the transformations q± (q± (q')) can lead only to one of the 
three momenta q', q+ (<?'), and q_ (q'). This property will be used below (see Tab.[T|. Since they lie outside of the real 
axis for q' = ±p n , the singularities of Z (q', q) play a role only in the case when q' corresponds to three free atoms. In 
the source term Z (q', qo) of Eq. (|45[) the singularities can appear only when both q' and qo correspond to three free 
atoms. However, the present work does not deal with this case of free-free transitions. Unlike the 3D case, where a 
partial wave expansion leads to logarithmic singularities of Z (q' , q), in the present ID case they are simple poles. 

For further derivation let us introduce dimensionless variables. Convenient energy and momentum scales are based 
on g, a characteristic value of the channel coupling strengths g m . The dimensionless momentum and energy are 
expressed as 

C=^(mg)- 2/3 q, e = m-^g-^E. (70) 

An introduction of new functions in place of the transition amplitude, the free Green function, and the two-body T 
matrix, 



F(CCo) = ~ m 7 /V /3 «oVe + «> - C 2 * (ff, 9b) (71) 
S (C',C) = —rnVW'^We + iO-FZ^q) = - — j-^ (72) 



r (0 = (e + i0- CT 1/2 T iD (k (?)) = ~ (uj^/e + zO - C 2 + |) \ (73) 



-l 



allows to exclude the square-root singularity from X (q,qo)- Here 

C: 



(O=-U±^Ve + i0-C, U eS = m^g- 2 / 3 U eS . (74) 



2^ 2 
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TABLE I: Algebraic properties of transformations £± 





c- (c- (O) 


c+ (c- (O) 


C- (C+ (CO) 


C+ (C+ (CO) 




c 


c+ (co 


c 


C- (CO 


-y/e/2 < (' < 


c 


C+ (CO 


c- (co 


c 


C > 


C+ (C) 


c 


C- (CO 


c 



Equation (|45p is transformed then to the form 

F (CO Co) = S (CO Co) + / ^(COC)r(O^(CCo). (75) 



A. Singularities of the free Green function 

According to Eq. (|72j) S (CO Q has simple poles on the real axis whenever C' 2 < £• The pole contributions can be 
separated using the symbolic expression 

' v.p.— ^-+i7n5(C-Co), (76) 



C - Co ~ M C - Co 

leading to the integral equation 

oo 

F (CO Co)-^r (C+) F (C+, Co)-i*r (C-) F (C-, Co) = -^—--^—+ [ d( v.p. f-i— - ) r (C) F (C, Co) • 

Co-C+ Co — C- j VC-C+ C-C-/ 

— oo 

(77) 

Here C± = C± (CO, symbol v.p. means Cauchy principal value integral in the vicinity of C±, and the source term is 
non-singular for the problems considered here, as it was mentioned above. The following derivation is based on the 
algebraic properties of the transformation C± (0- A direct calculation of products of these transformations (see Tab. 
H]) demonstrates that C± (C±) leads again to one of the three values CO C+, or C-- Therefore, replacement C' by C+ or 
by C- in Eq. (|77|) leads to two additional linear equations in F (CO Co), F (C+, Co), and F (C-, Co)- Then they form a 
system of linear equations with the matrix 

(1 — iirr (C,+) — «7rr(C-)\ 
-i7ir(C0 1 -wrr(C-) . (78) 
— iirr(C) ~* 7rr (C+) 1 / 

Solution of this system leads to the integral equation in F (CO Co) 

oo 

F(COCo) = S (CO Co) + J d(v.p.S(£,Or(()F((Xo). (79) 

— oo 

Here 

and the parameters s, s+, and s_ are expressed as 

8 = (A- 1 )^ sign (C + v^/2) + (A- 1 )^ sign ((> - y/Z/2) 

s+ = (A- 1 ) n - (A- 1 ) 13 sign (C - Vi/2) (81) 
S - = -(A- 1 ) 11 -(A- 1 ) 12 sign(C' + v ^/2) 



in terms of the inverted matrix (1781) . 



12 



B. Singularities of two-body T matrix 

Following derivation is a generalization of the approach [24| to the case of multiple bound states. Since r (() is an 
even function of £ [see Eq. (|T3|)] and S {—(', C) = S ((' , — C) (this property can be directly proven), the equations for 
the odd and even components, 

F±(C,Co) = F(C,Co)±F(-C,Co), (82) 

can be separated as 

oo 

F± (C, Co) = S ± (C, Co) + Jd( v.p.5± (C, C) r (C) F ± (C, Co) , (83) 

o 

where 

s ± (C,0 = s(C,0±s(-C,0- (84) 

The odd and even component in ID scattering are analogs of 3D partial waves [55[ . 

Equations ([83]) still contain singularities related to poles of r (C) at = ( n + iO, where Cn = ^ (wg)~ 2 ^ 3 p n . In the 
pole vicinity the function r (C) has the limiting behavior 

r ® r % 7 -n ' r " = * ReS ^ (0 - C - Cn) = ~ 4/ f 4 /3 -fi ( 85 ) 

C^Cn C - Cn - *0 ^ Tfl ' Q ' Cn 

[see Eqs. (|B2[) and (|73|) ]. The contributions of the singularities can be separated using Eq. (176)) leading to the 
regularized equations 

oo 

F± (C, Co) = S ± (C, Co) +nJ2r n S ± (£ , C„) F± ((„, Co) + v.p. / d(S ± (C, C) r (C) F± (C, Co) , (86) 

o 

where the Cauchy principal value integral is taken for all singularities in S± (C', C) an( i r (0- 

Source terms in these equations have forms of linear combinations of the functions S± (C')Cn) (f° r the problems 
under consideration Co is equal to one of Cn)- Thus the solution can be expressed as a linear combination 

F± (C, Co) = F ± (C, Co) + tt^ r n F± (C nj Co) F± (C, Cn) (87) 

n 

of solutions F± ((", Co) of integral equations 

oo 

F± (C, Cn) = S ± (C, Cn) + v.p. J d(S ± (C, C) r (C) F± (C, Cn ) , (88) 



where the source terms do not contain unknown functions. Solving of these equations also allows the evaluation of 
the values of F± at specific points, F± (( n ,(o), involved into Eq. ([57)) . They are determined by the system of linear 
algebraic equations 

Ft (Cn, Co) ~ 7T T n'F± (Cn, Cn' ) Ft , Co) = F± (Cn, Co) , (89) 

n' 

which are obtained by setting £' = Cn in Eqs. ([57)1 . 



Conclusions 



Multichannel resonant interactions can be described by a model involving one atomic and several molecular states. 
Even in the case of zero-range interactions this model predicts several two-body bound states (diatoms), which are 
superpositions of two-atomic and molecular components. Orthogonality of the diatom wavefunctions is proven. 
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In three-body problems the multichannel behavior leads both to an effective energy-dependent interactions and 
multichannel asymptotic wavefunctions. The last effect leads to rescaling of the probabilities of elastic and inelastic 
atom-diatom collisions and of three-atom association. All the probabilities can be expressed in terms of the symmetric 
transition amplitude for an effective single-channel problem. This amplitude can be evaluated as a solution of the 
Faddeev-Lovelace equation. Singularities in this equation can be regularized using their algebraic properties. The 
regularization procedure is applicable to generic ID three-body problems with separable interactions. 

The algorithm of solution of a resonant ID three-body problem consists then of the following steps: 

— functions F± (0 Cn) are determined as solutions of the regularized integral equations for each bound state 

n; 

— the system of linear equations (JHSJ) is solved for F± (Cn, Co); 

— F± (C, Co) is calculated by using of Eq. ([57]): 

— X (q,qo) is determined using Eqs. l[7Tj) and (|52"j) : 

— transition probabilities are evaluated using Eqs. (|49[) ([56]) . (|66|) and (|68|) . 



APPENDIX A: ORTHOGONALITY OF THE DIATOM STATES 



Integration of Eqs. (O and (|10p over an infinitesimal interval including z — allows to represent them in the forms 

,("0 (Al) 



1^(°) = eW°) (0) + V2]T 5m ^ 

ra 

±-d<pW=U eS (E n )ipW (0), 



(A2) 



(A3) 



where the wavefunction derivative discontinuity operator d acts as 

* dtp dip 

Clf — -^-2=0+0 7^-\z=0-0- 

oz oz 

Equations (|A1[) and ([6]) allow to prove the identity relation involving wavefunctions of two diatom states, n and n' , 

(0) ^ ] + E n £ ^ m) Vl m) = ^ (0) ^> + E n , ^ (A4) 

m m 

Multiplication of Eq. (|10| for tpi ^ by (pffl* with following subtraction of the complex conjugate of the similar product 
with exchanged n and n' leads to the following expression 



{E n -E n ,)^>{z)^{ Z ) = -l§- z 







J0) 



(0) 



9 (0)* 



d 



dz 



+5 (z) [U eS (E n ) - U cS (E n ,)] V W* (0) ^ (0) 
Integrating this expression from — oo to oo and using Eqs. (|A4[) and (|A2|) one gets the equation 



0. 



(A5) 



(A6) 



It leads to the orthogonality conditions ((15)1 for non-degenerate states, while the degenerate states can be always 
orthogonalized by a simple linear transformation. 



APPENDIX B: RESIDUES OF THE TWO-BODY T MATRIX 

For scattering on a generic potential U, the two-body T matrix, as a function of the collision energy E, has poles 
at the bound state energies E n with the residues 

Res(T 1D ,E = E n ) = U\^ n )(vn\U. (Bl) 

In the present case U = U e g (E) S (z) and the bound state wavefunction is given by Eq. (fTT)) . Taking into account 
Eqs. ([II]) and (jTSj) one obtains 

Res (Tin (k),k = w n ) = -2i-^W„. (B2) 

m 
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APPENDIX C: RELATIONS BETWEEN TRANSITION AMPLITUDES IN THE TWO-CHANNEL AND 

DEACTIVATION CASES 



The three-body asymptotic functions xidp d (qi, ki), corresponding to the deactivation product channels, are obtained 
by substitution of Eqs. ([M]) and (|3l)|) into Eq. (|42[) . For large \D&\ they are approximately independent of the 
momentum fe;, 



Xid Pd (qi,h) 



1/2 



6 (ffl - Pd) 



(CI) 



(this means that sizes of the deactivation-product bound states are negligibly small compared to other relevant scales) . 
The transition amplitude (fM)) can be then transformed, using Eqs. (f2l))) and (|4"3"|) , to the form 



I' dpd,lnp n 



J_ ( 2Wd 

2n \ 7TK d 



1/2 



dk v k / ^C n (ft",V) + 2(|) 1/ %^£ n (^) 



(C2) 



where n = 0, 1 and are solutions of the two-channel problem, which does not involve deactivation product states. 

Changing the integration variable ki> by qyi allows us to recognize in the integrals the coordinate-representations of 
the open and closed channel wavefunctions [see Eq. (|20p ] in the origin. Finally, the symmetric transition amplitude 
(|46| involving the deactivation product states can be expressed using Eq. (|44| as 



X (±p d ,p n ) = 3 



1/2 3 
K 



U*^ (0,0,0) + \/ -.MiM 



(C3) 



in terms of the values at the origin of the symmetric wavefunctions in the two-channel case 

^„(0,0,0) = 3- 1 / 2 ^vS n (0,0,0)= ^ (2n)- 3/2 J2 J dq 3 dk^ lnPn (q^kj) 



(C4) 



(0, 0) = 3-V> £ ^ (0, 0) = £ E / 



do,- 



£-M0-3g?/(4m) 



dkjip[ npn (qj,kj) 



[see Eqs. (J20J), flU, and lp|]. 

These values can be also expressed in terms of the symmetric transition amplitude for the two-channel case in the 
following way. Let us write out the Lippmann-Schwinger equation for the effective single-channel problem 



(qj, kj) = x lnp ( qj ,k 3 ) 5 jt + (e + iO - ff 00 - Uf ) 1 £ Uf [i>ln V 



Here the Green function for the free j th atom and diatom is the integral operator 



(£ +J 0-i7 00 -l> cff 



-1 r _• 



9o(qj,kj) / dk'. 



6(k' 3 -k 3 ) + ±-T 1D (k( qj ))g a (^fcj) 



expressed in terms of the two-body T matrix (|13[) and the Green function for three free atoms, 

3o (q, k) 



E + l Q-^-^-k 2 
4m m 



(C5) 



(C6) 



(C7) 



The incident wave xinp (q, k) [see Eq. (|4"l2"|) ] is proportional to go (q, k) (it is a general property of separable potentials, 
see [HI, HU). Substitution of -0 given by Eq. (|C5[) into Eq. (|C4[) allows to express the integrals in Eq. (|C4[) in terms 
of the transition amplitudes for the effective single-channel problem. Some algebra leads to the values of symmetric 
wavefunctions 
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